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The fluctuations of a received radio signal due to fading are assumed 
to behave like the envelope of narrow-band Gaussian noise. Estimates 
of the distribution of the fade lengths for various depths of fades are given, 
and relations which may be useful in analyzing fading data are derived. 
A similar problem involving the separation of the intercepts of the noise 
current itself, instead of its envelope, is also discussed. 

I. INTRODUCTION 

When a radio transmitter sends out a steady sine wave of very high 
frequency the signal received far beyond the horizon is no longer steady. 
The envelope of the received signal fluctuates in an irregular manner 
and may go through several maxima and minima in a second. These 
fluctuations have been studied by many people.* For observations 
extending over fifteen minutes or so, the variations in the size of the 
envelope follow, approximately, a Rayleigh distribution. For observa- 
tions extending over a month, the distribution of the hourly medians 
is roughly normal in decibels with a standard deviation of 8 or 10 db. 
There are also seasonal and diurnal variations. 



* In this paper we shall make use of the data given by K. Bullington, W. J. 
Inkster and A. L. Durkee. 1 
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Here we shall be concerned mostly with the short-term Rayleigh 
fading. In particular we shall be interested in the distribution of the 
lengths of fades, i.e;, in the lengths of the intervals during which the 
signal remains below some specified level — the lower the level, the 
"deeper" the fade. 

We shall assume, in our calculations, that the envelope of the re- 
ceived wave is the same as the envelope of a Gaussian random noise 
current having a narrow, normal-law power spectrum centered on the 
transmitter frequency.* Under suitable conditions the distribution of 
the fade lengths approaches a limiting form for deep fades and the exact 
form of the power spectrum becomes unimportant. 

Although the results obtained from the Gaussian noise model do not 
always agree with the rather meager available experimental data, they 
do show general trends. Moreover, a process of averaging over the 
parameters given in Section IV suggests relations between the depths 
of fades and their durations which may be of use in analyzing fading 
data. 

Our method of studying the duration of fades may be described as 
follows. We first find an expression which approximates the distribution 
function in the region of short intervals. This is done by using the method 
developed in Section 3.4 of Ref. 2 for dealing with the intervals between 
the zeros of a random noise current. The second step is concerned with 
the very long intervals. Although it seems to be generally accepted 
that the probability of a fade exceeding a great length T decreases as 
exp (— aT), no convenient method of calculating the exact value of a 
is available as yet. Here we more or less ignore this difficulty by using 
a method which amounts to plotting, on suitable coordinates, our ap- 
proximation for short intervals and then extrapolating by eye so as 
to obtain a "reasonable" looking curve. This method is described in 
Section VI. It turns out that the slope of the extrapolated curve at 
the origin gives the value of a. Needless to say, the reader should regard 
these portions (the long fade portions) of our distribution curves merely 
as estimates. 

In Section II we state some known results regarding the frequency 
and average duration of fades based on the Gaussian noise model. Sec- 
tion III lists the main results obtained in the later sections. In Section 
IV some of the results set forth in Section III are applied to radio 
fading problems and comparisons with experimental data are given. 
Section V is concerned with the lengths of the intervals during which a 

* This power spectrum is obtained when the propagation is assumed to result 
from randomly moving scatterers having a Maxwcllian distribution of velocities. 
Even if the power spectrum is not normal most of our formulas will still hold, 
but some of our curves will no longer apply. 
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Gaussian noise current I(t) exceeds some given value /. Although the 
material in Section V is not directly related to the fading of the envelope, 
it is of interest in its own right and serves to set the stage for the work 
of Section VII, which deals with the intervals during which the envelope 
R(t) of l(t) lies below some given value R. A discussion of the computa- 
tions is given in Section VIII. Several known results concerning the 
.separation of the zeros of a simple Markoff process are listed in Section 
IX and their relation to the present work is pointed out. 

II. THE AVERAGE DURATION OF A FADE 

Here we shall state several known results regarding the fluctuation of 
a Gaussian noise current /(/) and its envelope R(t)* Let the mean square 
value of I(t) be b (b is the same as ^o and b of Refs. 2 and 3). Then 

probability 1(1) > I = |[1 - P(Ib~ m )], (1) 

probability R(l) > R = e"' m ", (2) 

where P(x) is the error integral 

™-^f «-""«<• (3) 

Let w(f) be the power spectrum of 1(f) and /3 be the — \pa" of Ref. 
2 so that 

£ = f w(f)C2-KJf df. (4) 

Jo 

Then, according to a result due to Mark Kac, the expected number of 
times per second I(t) passes upward across the level I(t) = I is 

Nl = i_ G3/6) l V J/26 . (5) 

Air 

The average length of the intervals during which /(/) > / is the quo- 
tient of (1) and (5): 

t = ir(b/ii) ll2 e' ll2b [l - P(Ib~ m )\, 

E~|(2ir//3) 1 ' 2 , I»b m (G) 

t~ 2t(&//8) 1 'V , '* > I«-b m 



* See Refs. 2 and 3, where reference is made to the earlier work of V. D. Landon, 
K. A. Norton and others. 



584 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1958 

where the first asymptotic value follows from 

Ml"Hr)]~^4=' *»1. (7) 

The corresponding results for the envelope run as follows.* Let w(f) 
be confined to a narrow band and be symmetrical about the midband 
frequency / , and let /3 r be the b 2 of Ref . 3 so that 

0, = 4/ f w(fKf - /o) 2 df. (8) 

Jo 

Then the expected number of times per second R(t) passes downward 
across the level R is 

N R = ^ r /2^'\R/b)e- RV2b . (9) 

The average length of the intervals during which R(t) < R is, from 
(2) and (9), 

t = (b/R) (2TM il2 (e R:,2b - 1). (10) 

When R becomes small 

t -> (fi/2) (27r//3 r ) 1/2 . (11) 

In most of the following work we shall set b equal to 1 so that 

pffi = b = f w(f) df = 1. (12) 

Jo 

When this is done we have 

t/t = [1 - P(I)] exp (7 2 /2), 

tf,/tfo = exp (-/ 2 /2), 

where N denotes the average number of upward crossings, per second, 
across the axis I(t) = and t the average separation of the zeros of 
I(t). Table I gives values of these ratios. They do not depend upon the 
power spectrum w(f). 

The ratios corresponding to (13) for R(t) are 

tf tl = (e™ - 1)/R(e m - 1) « 1.54 (e™ - 1)/R, 

(14) 
Nn/Ni = R exp [(1 - R 2 )/2] « 1.649 Re R ''\ 



* K. A. Norton and his colleagues 4 have considered more general fading prob- 
lems in which the received wave contains both a steady sinusoidal component and 
a Rayleigh-distributed component. 
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Table I — Average Durations and 
Crossing Rates — I(t) 
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7/rms /(/) 


t/k 


Nl/Nt 


Nft 


-5 


536,674 


3.73 X 10"» 


1 - 2.87 X 10- 7 


-4 


5,962 


3.35 X 10~ 4 


0.999968 


-3 


180 


0.0111 


0.9986 


-2 


14.4 


0.135 


0.977 


-1 


2.77 


0.607 


0.841 





1.000 


1.000 


0.500 


1 


0.523 


0.607 


0.159 


2 


0.336 


0.135 


0.0228 


3 


0.243 


0.0111 


0.00135 


4 


0.189 


3.35 X 10" 4 


3.17 X 10- 6 


5 


0.154 


3.73 X 10-« 


2.87 X 10~ 7 



where ~t\ and iVi refer to the level 72=1. Other reference levels some- 
times used are 



median R = (2 log 2) 1/2 = 1.1774 
average R = (x/2) 1/2 pa 1.25, 
rms R = 2 m Pa 1.41. 



1.18, 



(15) 



Thus, instead of the second of equations (14), it is sometimes convenient 
to use 

R-/R m - 



Nr = 2R _« 2/2 = 2R/l\ 
N m R m e R m \2/ 



where N m and R m refer to the median level. In (14) and (15) the value 
of b is unity, i.e., rms I(t) = 1, I(t) being the noise current whose en- 
velope is R(t). Table II gives values of the ratios (14). Just as in the 
case of I(t), these ratios do not depend on the power spectrum. 

Table II — Average Durations and 
Crossing Rates — R(t) 



R/ rms 7(0 


<7<t 


Nr/Ni 


NRt 


0.125 


0.0967 


0.2045 


0.00778 


0.25 


0.1957 


0.3995 


0.0308 


0.5 


0.4105 


0.7275 


0.1175 


1.0 


1.000 


1.000 


0.3935 


1.18 (median) 


1.309 


0.9731 


0.5000 


2.0 


4.924 


0.4462 


0.8647 


3.0 


45.74 


0.0550 


0.9889 


4.0 


1148 


0.002212 


0.999665 
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III. GAUSSIAN NOISE RESULTS — NORMAL-LAW POWER SPECTRUM 

Let 1(f) be a Gaussian noise current with the power spectrum 

„,(/)= —7T<r' w (16) 

<T V27T 

so that the mean square value of I(t) is b = 1. When I(t) crosses upward 
through the level I the probability that it will remain above I for more 
than t seconds may be denoted by F(u, I). Here 

U = r/t, 

t = (2a)- 1 [1 - P(I)\ exp (/ 2 /2), (17) 

P = W, 

where t is the average length obtained from (6). 

Fig. 1 shows our estimates of F(u, I) plotted against u for various 
values of 7. These curves were obtained by the calculations described 
in Section VIII. As mentioned in the Introduction, the values of F(u, I) 
for the larger values of u must be regarded as conjectures. Two special 
cases are 



F(u, /-> * ) = e 



< I8 > 



and these hold not only for (16) but apparently for any power spectrum 
such that F(u, I) exists. The gradual change of F(u, I) from an exponen- 
tial form to a Rayleigh form as / runs from -co to + °° [as indicated 
by (18)] was surmised by D. S. Palmer from an examination of some 
experimental data.* 

A case which has received some attention in the literature occurs 
when 

w(f) = 8/[l + (27T/) 2 ] 2 . (19) 

This power spectrum is unusual because it represents a borderline case 
between a simple Markoff process, for which F(u, I) does not exist, and 
the smoother noise processes, one of which is represented by the w(f) of 
(16). An experimental investigation of the spacings between zeros of 
I(t) when w(f) has the form (19) [and also when w(j) has several similar 
forms] has been carried out by R. R. Favreau, H. Low and I. Pfeffer. 6 

* Mr. Palmer has informed me that the "Note added in proof" appearing in 
his paper should he disregarded. He has ohtained a proof that the Rayleigh form 
is indeed approached at I = °° . His proof is given in our Section V under the 
discussion of equation (66). 
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Fig. 1 — The probability F(u, I) that 7(0 > / for an interval lasting longer 
than t. /(/) has a Gaussian power spectrum and I is given by equation (17). 



D. S. Palmer 5 has obtained an approximate integral equation for the 
probability density /(«, 0) = — dF(u, 0)/du which has been improved 
by J. A. McFadden. 7 McFadden gives curves obtained by solving 
his equation numerically. They agree well with the experimental re- 
sults of Ref. 6. For (19), /(w, 0) appears to approach a value close to 1.25 
as //approaches zero, whereas for (16), /(m, 0) approaches zero. 

The corresponding results for the envelope R (/) may be stated as fol- 
lows. Let /(/) be a Gaussian noise current with the envelope R(t) and 
the narrow-band power spectrum 



w(f) = <T\2*r m exp [-(/ - /o)72a 2 ]. 



(20) 
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When R(t) crosses downward through the level R, the probability that 
it will remain below R for more than t seconds may be denoted by 
F T (u, R), where now 

u = r/t, P r = (2ttct) 2 , 

t = (Ra)- 1 (2i r )- m (e Rm - 1), 



(21) 



and I is computed from (10). 

Estimates of F T (u, R) are given in Table III and the calculation of 
these values is discussed in Section VIII. Just as for 7(0, the values for 
the larger w's must be regarded as conjectures. The special cases corre- 
sponding to (18) are 

(22) 



F r (u, R-^0) = (2/u)I 1 (2/iru 2 ) exp (-2/ttw 2 ), 



(23) 



where Ii(z) denotes a Bessel function of imaginary argument. Again these 
limiting forms do not depend on w(f). The values given in Table III are 
plotted in Fig. 2. It may be shown from (23) that, for small values of 
u, F T (u, R — > 0) is approximately 1 - ^wu and, for large values of u, 
it decreases like 2/tm\ Since the t in the u of (23) is given by (11), u is 
proportional to t/R, and hence F,(u, R — » 0) depends on t and R only 



Table III — Estimates of F,(u, R) For the Narrow-Band, Normal- 
Law Power Spectrum (20) 



« 


R = 


R = 0.5 


R = 1 


R = 1.18 


R = 2 


R = oo 





1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.1 


0.995 


0.994 


0.993 


0.992 


. 0.976 


0.905 


0.2 


0.976 


0.974 


0.969 


0.963 


0.885 


0.819 


0.3 


0.945 


0.938 


0.924 


0.906 


0.793 


0.741 


0.4 


0.894 


0.884 


0.843 


0.80S 


0.711 


0.670 


0.5 


0.823 


0.804 


0.733 


0.694 


0.637 


0.606 


0.6 


0.729 


0.697 


0.624 


0.598 


0.571 


0.549 


0.7 


0.619 


0.588 


0.535 


0.528 


0.511 


0.496 


0.8 


0.521 


0.492 


0.464 


0.467 


0.458 


0.449 


0.9 


0.429 


0.412 


0.407 


0.413 


0.411 


0.406 


1.0 


0.354 


0.343 


0.360 


0.364 


0.368 


0.368 


1.2 


0.242 


0.249 


0.279 


0.285 


0.296 


0.301 


1.4 


0.171 


0.190 


0.216 


0.222 


0.237 


0.246 


1.6 


0.122 


0.149 


0.168 


0.174 


0.190 


0.202 


1.8 


0.0890 


0.116 


0.130 


0.135 


0.153 


0.165 


2.0 


0.0681 


0.0899 


0.101 


0.106 


0.123 


0.135 


2.2 


0.0535 


0.0701 


0.0784 


0.0826 


0.0987 


0.111 


2.4 


0.0425 


0.0544 


0.0608 


0.0645 


0.0792 


0.0907 


2.6 


0.0340 


0.0421 


0.0472 


0.0504 


0.0636 


0.0743 


2.8 


0.0269 


0.0327 


0.0366 


0.0393 


0.0510 


0.0608 


3.0 


0.0220 


0.0253 


0.0284 


0.0307 


0.0410 


0.0498 
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1.2 1.4 1.6 
U=7/t 



Fig. 2 — The probability F r (u, R) that R(t) < R for an interval lasting longer 
than r. R(l) is the envelope of narrow-band Gaussian noise having the power 
spectrum (20). 

through the ratio t/R. The curves for F r (u, R — > 0) and F r (u, R — ► oo) 
cross near u = 5.68, where they both have the value 0.00340. This be- 
havior is shown in Fig. 4. 



IV. COMPARISON OF THEORY WITH EXPERIMENT 

By combining the average fade lengths given by Table II and the fade 
length distributions F r (u, R) given by Table III, we obtain Fig. 3. In 
Fig. 3 the scales have been chosen so that the median value of R(t) is 
1.18 [see equations (15)], and so that fades below the median last one 
second, on the average. This leads to t x = 0.763 seconds and /3 r = 4.53 
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Fig. 3 — Distribution of the duration of fades, ltd) is the envelope of narrow- 
band noise. The time scale is chosen so that fades below the median value of R(t) 
have an average length of one second. 

and makes a = 0.339 in the associated power spectrum (20). When both 
R and r/R are small, the ordinates of the curves on Fig. 3 are given ap- 
proximately by 1 — 1.7(r/B) 2 . As R increases, the distribution changes 
from the form shown by (23) to the exponential (22). 

When the available experimental results for rapid fading are studied, 
it is found that some of them seem to agree with Fig. 3, but others do 
not. It is not hard to find reasons for lack of agreement when it does oc- 
cur. One is that we have arbitrarily taken the normal-law power spec- 
trum (20) for our calculations. This gives fade length distributions 
F r {u, R) of the sort shown in Fig. 2. These distributions appear to be 
typical of those obtained for w(f)'s which (i) decrease rapidly for fre- 
quencies outside the nominal band and (ii) do not become excessively 
large for frequencies near/o (such frequencies correspond to slow fluctua- 
tions in the envelope). For w(f)'s which do not satisfy these requirements, 
say one which behaves like 



w(f) = 



2/,- 



(24) 



A-' + -±7T 2 (/ - /o) 2 

when / is near / , k being very small, the work in Appendix IV suggests 
that the distribution curves have the quite different shape shown in Fig. 
13. When the experimental curves differ from those of Fig. 2, they tend 
to look like those of Fig. 13. 
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Fig. 4 — Plots of F r (u, R) for special cases. F r (u, R — > =■= ) is given by (22). 
F r (u, R — » 0) corresponds to fixed values of the median and fading rate, and is 
given by (23). 4>(r/i?) of (30) corresponds, to fixed median but variable fading 
rate . 



A second reason for the lack of agreement is that Fig. 3 refers to 



fixed values of the parameters b and f3 r where b ' is a measure of the 
median value (R m = 1.186 ") and (/3 r /6) " is a measure of the rate of 
fluctuation [the average number of times R(t) passes down across the 
median value in one second is N , fm = 0.235 (/3 r /&)' 2 ]. Actually, assuming 
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the Gaussian noise model to be valid, the values of & and /3 r may change 
slowly during an experimental run. When the fluctuation is regarded as 
due to the Doppler effect of moving scatterers, (/3 r /6) 1/2 is proportional 
to the average speed of the scatterers, which we might expect to be 
roughly proportional to the wind velocity in the scattering region. The 
value of b is the average received power and is proportional to the num- 
ber of scatterers and to the average power scattered to the receiver by 
each. According to this picture, it is not surprising that /3 r and b may 
change during an experimental run. 

In spite of all of these uncertainties we shall attempt (perhaps rashly) 
to draw some conclusions from our formulas in the case of deep fades 
(R near zero) . In this case a few of our formulas take a relatively simple 
form, and we may average over values of b and p r or otherwise use the 
formulas as a guide. The reader should bear in mind that the statements 
in the remainder of this section are tentative and future experimentation 
may show them to be in error. 

As an example of the averaging process, let us consider the expected 
number of times per second, N R , the envelope of the received wave 
crosses downward across the value of R microvolts where R is very small. 
When we assume b and /3 r to have the joint probability density g(b, /3 r ) 
and take the average of expression (9) (with R so small that the ex- 
ponential term is almost unity), we find that if g(b, /3 r ) satisfies suitable 
conditions, 

N R ttR f J V l {fi r /2v) m g(b, ft) db dfr (25) 

= R&. 

Here Ci is a constant which depends on the probability density of b and 
/3 r . By dealing with the probability that R(t) < R in much the same way 
and then dividing by (25), we can show that the average length of the 
intervals during which R (t) < R is 

t « RC 2 (26) 

where d is a constant of the same nature as d and R is again small. 

Equation (25) says that, for a given experimental run, the average 
number of times per minute the signal drops below R/2 microvolts 
should be half the number of times it drops below R microvolts, R being 
very small. Furthermore, from (26), fades below R/2 microvolts should 
last, on the average, only half as long as the fades below R. These rela- 
tions should hold for small values of R even though the experimentally 
observed ratios for the larger values of R may not agree with Table II. 
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The same sort of argument may be applied to the fade length distribu- 
tion. As R becomes small, F T (u, R) approaches a definite function of 
u = T /t, say <p(ii), provided u is not too large. Let R have the same sig- 
nificance as in (25). Then the probability that R(t) will remain less than 
R for an interval lasting longer than t seconds is approximately 

& I C r " 4/I {1 1/I] »»■ M db * - * (r/B) - (27) 

where iV* is given by (25). Thus the probability depends on t and R 
only through the ratio t/R when 72 is small. The curves of Fig. 3 corre- 
sponding to the smaller values of R show this behavior (as indeed they 
should). For example, the curve for R = 0.125 is obtained almost exactly 
by shifting the curve for R = 0.50 to the left by an amount correspond- 
ing to quartering the fade lengths. 

As an illustration of how the averaging process (27) may change 
F T (u, R), we suppose that the median level remains fixed at R = 1.18 
while the rate of fluctuation changes in such a way as to formally carry 
the normal-law power spectrum (20) over into the power spectrum. 

w(f) = Trfc 2 [h 2 + 4tt 2 (/ - /o)T 3 ' 2 . 

This particular example is chosen so as to simplify the steps which fol- 
low. Instead of introducing g(b, r ) and averaging over /3 r it is more 
convenient to average over the a of (20). Some mathematical experi- 
mentation indicates that we should take 

(/c 2 /4irV) exp [-fc78*V]der (28) 

to be the probability that a lies in the small interval da because this 
choice makes the integral from a = to a = =o of the product of (28) 
and (20) give the required w(f). When we give r the value (2ir<r) 2 in- 
dicated by (21), the <p function in (27) becomes 

*V!E] (29) 

where we take <p(u) to be the limiting form of F r (u, R — > 0) given by (23). 
The integral (25) now gives N R = Rk/2 and the integral in (27) becomes 
the integral from <r = to <r = <*> of the product of (28) and (29). The 
integration can be performed and gives 

Hr/R) = [l + '{]"" [l + ^ + « |/l + f ]"', (30) 

where u = kr/R. This is the fade length distribution, when R is small, 
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which results when the median level remains fixed but the rate of fading 
fluctuates according to the law (28). The average value t of the fade 
length works out to be R/k, so that u = t/1 in (30), in agreement with 
our earlier notation. 

If we were to take a" in place of a in (28) (and change the multiplying 
constant) we could obtain a more general, but also more complicated, 
expression for $>(t/R). When v — 2 the power spectrum (24) is obtained 
but the integrals (25) and (27) corresponding to N K and $(t/R) do not 
converge. 

The expression (30) is plotted as a function of u in Fig. 4, along with 
Friii, R —> <*>) and F r {u, R — > 0), given respectively by (22) and (23), 
for comparison. As u increases, (30) ultimately decreases as 2/w 8 . It is 
seen that the $(t/R) given by (30) and F r (u, R — > 0) [from which (30) 
is derived] arc quite different in appearance for small values of u. The 
distribution for fixed median and fixed fading rate is F r (u, R — ► 0), 
while the $(t/R) of (30) corresponds to a fixed median but a variable fad- 
ing rate. It appears from (27) that many different &(t/R)'s [of which 
(30) is only one] may be obtained by averaging over the median values 
and fading rates with different weighting functions g(b, fi r ). At this stage 
we are unable to tell from either theory or experiment whether a dif- 
ferent $(t/R) will be required for each set of data, or whether there are 
only two or three universal types of 3>(t/JR)'s which will fit, for small R, 
all kinds of fading data. 

Suppose that we are given a set of experimental results on short-term 
fading. How can we go about analyzing them to see whether they have 
anything in common with our theory? 

One of the easiest things to test is how N R , the average number of 
fades per second below the level R, depends on R. According to (25), 
Nn should be proportional to R when R is small, and, if the short-term 
median value and the rate of fading have remained fixed throughout the 
experiment, the ratios of Nn/Ni should agree with those of Table II. 
From (20) the same is true of t and l/ti . However, when the short- 
term median value and the rate of fading fluctuate, N R and t should 
still be proportional to R when R is small, but the ratios in Table II may 
no longer hold. It should be noted that Table II lists values of both 
N „/ Ni and l/ti , and it may turn out that one set of ratios will agree 
better with experiment than the other. This may happen because N R 
and I do not depend on b and /3 r in the same way. N R depends on both 
/3 r and b, t only on /3 r , and the probability that R(t) < R (which, for 
small R, becomes R /2b) depends only on b. In fact, the deviations from 
Table II and from the Rayleigh distribution should give clues to the 
variation of b and (3 r . 
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When experimental fade-length distributions for several levels of R 
are available, it is instructive to plot them on semi-log paper as a func- 
tion of t/R in much the same way as F r (it, R) is plotted as a function of 
u in Fig. 2. As R becomes small the curves should tend to coincide, at 
least when t/R is not too large. This follows from (27). If the limiting 
curve can be made to coincide with the R = curve in Fig. 2 by chang- 
ing the scale of t/R, then it is likely that the other curves in Fig. 2 and 
the curves of Fig. 3 will be valid. This indicates that b and /3 r did not 
change much during the experiment. 

If, on the other hand, the limiting curve looks like those of Fig. 13 or 
like the $(t/R) of Fig. 4, the probability of long fades is much greater. 
In this case one should plot on semi-log paper R times the probability 
that the fade length exceeds t versus t/R 2 . For large t and small R, the 
points should tend to lie on a common curve, which becomes a straight 
line as t/T 2 becomes large. This behavior is suggested by one of the con- 
jectures of Appendix IV concerning long fades when w(f) behaves like 
(24) near / , k being very small (i.e., much power in the frequencies 
corresponding to slow fluctuations). It is assumed that w(f) is still such 
that t exists. For w(f)'s of this sort the conjecture says that when t is 
very large and R is small 

P(t, R) = Fr(it, R) -^ Ul-RT 2 E e~ ai " r (31) 

where P(t, R) is the probability that a fade below level R will last longer 
than t seconds. In (31) the median value of R(t) is 1.18 and on isj, 2 R~ 2 , 
ji being the ith zero of ./o (•?*). 

From (11), t is proportional to R and it follows that RP(t, R) tends 
to be a function of t/R 2 . This is the basis for the method of plotting 
described above. When r is very large, only the first term in (31) is im- 
portant, and it gives a straight line on the semi-log paper. As R increases, 
«! is no longer given by 5.78 R~ 2 . For R = 1, «i « 5.0 and for R = 
median = 1.18, ai « 3.25. 

Unfortunately we do not have any clear-cut rules to tell when, as t 
increases, the fade-length distribution changes from a function of t/R 
to \/R times a function of t/R 1 . 

Sometimes experimental results are given in the form of curves show- 
ing the average number of fades per hour, say H(t, R), below the level R 
which last longer than t seconds. According to (25) and (27), H(t, R) 
should take the form 3600 NMt/R) = CR$(r//?), where C is a con- 
stant, R is small, and t/R not too large. Under the conditions leading 
to (31), we expect H(t, R) to be of the form ^(t/R 2 ) when t is large 
and R is small. In fact, if we assume (31) and note that Nut approaches 
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R /2 [b is one in (31), so that the median value here is 1.18], 

H(t, R) -> 7200 k exp [-5.78 kr/R 2 ]. (32) 

Since this expression contains only one parameter, k, and is based on a 
number of assumptions, it is probably too simple to agree well with ex- 
periment. Nevertheless, if the experimental data for H(t, R) passes the 
test of being of the form ^(t/R 2 ) for large r and small R, it would be 
interesting to go one step farther and see how well the data can be 
fitted by a suitable choice of k in (32). 

It is interesting to examine Fig. 13 of Ref. 1 in the light of the above 
results. This figure, with some dashed lines added, is reproduced here as 
Fig. 5. The three solid lines give the observed (for 505 mc) number of 
fades per hour which were more than 5, 10 and 15 db below the hourly 
median level and lasted longer than the time intervals shown on the 
abscissa scale. The two dashed straight lines labeled (a) show the 5- 
and 15-db curves obtained by applying our results to the 10-db curve in 
the following way. The number of fades per hour of depth R of duration 
longer than r seconds is H(t, R) which, when R is small and t/R not 
too large, may be written as 

H(t, R) = CR*{t/R), (33) 

where C is a constant. The values of R corresponding to the 5-, 10-, 
15-db curves are in the ratio of 1.78, 1, 0.56. Since the curves of Fig. 5 
are plots of H(t, R) versus t, (33) gives a relation between them. Thus, 
suppose the coordinates (ti , Hi) of a point on the Ri curve are known. 
Then the coordinates of the point (n , H 2 ) corresponding to the same 
value of t/R on the Ri curve are given by 

T2 = TlRz/Rl } 

(34) 
Hi = H1R2/R1 . 

The 15-db dashed line (a) was obtained by taking the solid 10-db 
line to be the Ri curve and using (34) with R2/R1 = 0.56. The 5-db 
dashed line (a) was obtained in a similar way from the 10-db line by 
taking R2/R1 = 1.78, even though this may take us into the region where 
our approximation begins to fail. The dashed lines (a) are drawn only 
for the smaller values of t/R because of our assumption that (33) may 
not hold if t/R is too large. 

Under certain conditions we may expect H(t, R) to depend only on 
the ratio t/R 2 when t is large and R is small. The dashed straight lines 
labeled (b) are obtained from the 10-db curve by assuming H(t, R) to 
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Fig. 5 — Number vs. duration of fast fades at 505 mc. 

be of the form ^(t/R 2 ) and proceeding in much the same way as in the 
calculation of the lines (a). The lines (b) are drawn only for the larger 
values of t/R 2 because t is supposed to be large. 

The curved dashed lines labeled (c) in Fig. 5 are plots of H(t, R) 
computed from (32), with k chosen so that H(t, R) agrees with the 
experimental 10-db result at the ordinate corresponding to ten fades 
per hour. Since the value of R, 10 db below the median value of 1.18, 
is 0.372, and since t = 3.8 seconds corresponds to ten fades per hour, 
the equation to determine k is 

#(3.8, 0.372) = 10. (35) 

This yields two values of k, namely 0.0018 and 0.0147. The value 0.0147 
was used in plotting the curves (c) in Fig. 5, since the other value gives 
curves with much smaller slopes. 

The difference between the dashed curves in Fig. 5 and the correspond- 
ing solid curves is noticeable. An examination of the original data indi- 
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cates that some of the discrepancy may be due to the fact that, in one of 
the experimental runs, N K for JO db exceeds the N H for 5 db. This 
completely contradicts (25) and is probably caused by the short-term 
median (as measured by the parameter 6) temporarily dropping down to 
a low value. This violates our assumption of R being "small". The rela- 
tion between the 10-db and 15-db experimental results seems to agree 
fairly well with the theory, but comparison is difficult because of the 
scantiness of the 15-db data. 

Our results for R(t) may also be compared with some experimental 
results given by R. R. Favreau, H. Low and I. Pfeffer. Their work is 
concerned with the probability that in an interval of T seconds duration 
there will be at least one fade of depth R which lasts longer than r 
seconds. They express this probability as 1 — exp [— P'(0)f b T], where 
fb is the associated bandwidth, and their P'(0)fb may be roughly equated 
to our NrF t {t/I, R). Although the comparison is not straightforward, 
it appears that there is fairly good agreement for values of R larger 
than 0.25. As R decreases below 0.25 a disturbing discrepancy between 
experiment and theory seems to creep in. The reasons for this are not 
clear. 

Although in this paper we are primarily concerned with rapid fading 
and have been speaking about it in this section up until now, it is of 
some interest to try to apply our theory to the slow fades, i.e., to the 
distribution of the hourly medians. Since the hourly medians have a 
distribution which is nearly normal in decibels, it is natural to compare 
the observed statistics with those of a Gaussian noise current I(t). In 
this comparison "db above the median" plays the same role as "am- 
peres". 

When the fade-length distributions given in Figs. 6 and 7 of Ref. 1 
are plotted on semi-log paper the resulting points tend to lie on straight 
lines instead of on curves shaped like those of Fig. 1. This exponential 
distribution of fade lengths indicates that the power spectrum of I(t) 
is more like the ones discussed in Section IX and Appendix IV than 
like the normal type shown in equation (16). A stochastic process which 
does have an exponential distribution of fade lengths may be constructed 
by drawing a value at random from a normal universe every hour on 
the hour. However, the successive values in such a process would be un- 
correlated. This is not true for the hourly medians. 

Even though we have little justification for doing so, it is interesting 
to compare the separation between the experimental curves for various 
depths of fades with the separation predicted by the analogue of equa- 
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t.ion (27). The analogue states that the probability of 7(0 remaining 
less than — 7 for an interval lasting longer than t is approximately of 
the form $(tI) when / is large. According to this rule, the — 10-db 
curves in Figs. (> and 7 of Ref . 1 may be obtained by shifting the — 5-db 
curves to the left by an amount corresponding to changing the time 
scale in the ratio 5 to 10 (or 1 to 2). Actually, a ratio of 1 to 3 would be 
required to agree with the experimental results. The analogue of (31) 
suggests in a vague way that when t becomes large, / remaining large, 
the function $(tI) might change into a function of tI 2 . This would 
correspond to changing the time scale in the ratio of 1 to 4. However, 
not much significance can be attached to these figures. 

V. DISTRIBUTION OF THE LENGTHS OF THE INTERVALS DURING WHICH 
Kt) > I 

In Section II we have given expressions for the average length I and 
the expected number per second Ni of the intervals during which 7(0 
exceeds the value /. As a first step in finding the probability p(r, I) dr 
that the length of such an interval lies between r and t + dr, we obtain 
an approximation, pi(r, /) dr, which agrees well with p(r, I) dr for 
small values of r. The expression pi(r, 7) dr is the probability that, if 
7(0 passes upward through the value 7 at time t t , it will pass downward 
through / in the infinitesimal interval fe , h + dr where fe = h + t. 

An expression for p\(r, I) may be obtained by the method outlined 
in Section 3.4 of Ref. 2. Thus, the probability that I(t) will pass upwards 
through the level / in /i , /i + dt\ and downwards through / in t-> , (■> + dl-. 
is 

-dhdkj dli j_ dh'h'tipil, U, h', D, (36) 

where p(Ii , //, /•>' , !■>) denotes the probability density of 7(/i), and 
7'(/i), 7'(/ 2 ), 7(/ 2 ) and the primes denote time derivatives. Then 



, T\ i t = expression (36) 

' ' T ' " prob of 7(0 passing upward through 7 in U , h + dh 



(37' 



and we obtain 

Pl (r, 7) = -2wfir w f f dti I dh'h'h'pd, Ii, It, D- (38) 
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Here is given by equation (4) and we have taken b = I 2 (t) = 1. 
The probability density p(h , h', I 2 , I 2 ) is 

(2t)- 2 M- U2 exp I - ^ [MrJi + 2M n hU + 2il/ 13 / 1 / 2 ' 

+ 2M 14 /i/ 2 + Mnli* 4- 2M 23 h'l 2 



(39) 



+ 2M 2i I 1 'l 2 + M 3 3/ 2 /2 + 2Muh'h + M 44 / 2 2 ]j, 



where il/,y is the cof actor of the i/th element in the determinant 

10 mm 

r\ r> II I 

13 —m —m 

m —m /3 

m —m 1 

in which we have simplified the notation of Ref. 2 by writing 



M = 



(40) 



= ^r = / w(f) cos 2tt/t df, 
Jo 



(41) 
m! = dm/dr, m" = d 2 m/dr 2 , /3 = [— m"] T=0 . 

The values of the M,-/s and some relations among them are listed in 
Appendix I. When h and I 2 are set equal to I and use is made of the 
fact that some of the Af,-/s are equal, the exponent in (39) reduces to 
- 1/(2M) times 



2(M n + M l4 )I 2 + 2(M 12 + M u ){Ih' - lift + ¥„(//' + i/) 

+ 2il/ 23 7i7 2 / . 



(42) 



Replacing //, I 2 in (42) by Zi — ai , .r 2 — a 2 , choosing the o's so as to 
make the first order terms in Xi and x 2 vanish and using 



gives 



(M 12 + 3/24) (1 + m) = m'(M 22 - M n ) 
M = (M n + il/ 14 )(l + m) - m'(M 12 4- M M ) 

ai = — a 2 = m'I/(l 4- m) 



(43) 



(44) 



and converts (42) into 
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il/ 22 (.ri 2 + x?) + 2M^c 1 x i + 2I 2 M(l + m)-\ (45) 

When the change of variable corresponding to 

x(l - r 2 ) -1 ' 2 - Xi(Ma/M) w 

y{i _ r *)-W = -x 2 (M 22 /M) m (46) 

w (1 - m ) + mm 
r = Mn/M n = — - — 

0(1 — m 2 ) — m 

is made in (38) it goes into 

p l(Tj /) = Mn0- m (l - m 2 )' 3 ' 2 exp [£ - -iL] /fr.fc), (47) 

where 

./(,., /,) « JL ( <fc ( rfrfz - fe)(j/ - h)e% 

ZlTS Jh >>h 



_ m7_ \lzjlj\ s _ (1 - ,-=)'«, 
1 + m |_ Af 22 J 

M 22 = /3(1 - w 2 ) - m , 

m _ ,"- +jg - 2r,, = _ 2 _ V 2 (t 2 + ^2 _ 2?<x2/)> 
2(1 - r 2 ) 



(48) 



(49) 



and we have used 

(1 - m*)M = M n 2 - M a * = Mn\l - r"). (50) 

J(r, h) may be expressed in terms of the tabulated double integral 

K(r, h) = J- I dx [ dye' (51) 

Z7TS JA A 

as follows. We have 

6 z —2/ \ z 

dx 

d t —2/ \ i 

— e = s (rx - y)e . 
dy 

Elimination of the y terms on the right leads to 

(£+'£)'- -- r - (52) 
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In much the same way, elimination of the x~ and y" from the second 
partial derivatives gives 



r n + r Fl + C 1 + r '^ ^r 
ox* ay- oxay_ 



e~ = xye' — re". (53) 

When the integrand in J(r, h) is written as (xy — hx —hy -\-h)e z 
and use is made of (52), (53) and the symmetry in x and y, we obtain 

»2 



2rf,+ U + r 2 )-V 

ay- axay 



J(r, h) = J - f </.<■ f dy 

+ r + 2/K1 + r) i- + fe 2 ] e' (54) 

= (;• + /r)/v (/-,/,) + (a/2,) exp f- J^-A - fe(2 T )- 1 V fc ' /8 [l - P(a)] 



where P(a) is the error integral defined by (3) and 



•-<fe]" r 



/r 



a 2 + h 2 



+ r 



(55) 



Values of K(r, h) may be obtained from tables 8 of d/N. This is true 
because d/N is the function of the three variables r, h, k obtained from 
the right-hand side of (51) by changing the lower limit of integration 
in the y integral from h to k. The tables do not extend to negative values 
of h, and K(r, h) for these values must be obtained from 



K(r, -h) = K(r, h) + P(h). 
From (56) it may be shown that 



M, 



Pi(t, -I) = Pi(t, /) + J 



j gl/2( 1 _ OT 2)8/2 _ 



2Ae" fc,/ H 
V2ir J 

exp B - rrJ 



(56) 



(57) 



where h depends on I in accordance with (48). 

The following special values of K(r, h) and J(r, /i) may be derived 
directly from known results. 

1. When h = 



, v 1 cos ' r 
K(r, 0) - - - — - — , 



< cos r £ 



(58) 



J(r,0) - rK(r,0) + 



2tt' 
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2. When r - 



s = 1, a = h, 

K(0, h) = i[l - P(h)f, 
_ \h , __ „ /lA1 



(59) 
./((),/,) = Jjj[l -P(fc)]- (27r)- 1/ V 

3. When r approaches unity from below we have, in the limit, 

s = 0, o = 0, 
A'(l, h) = |[1 - P{h)] t (60) 

.7(1, h) = 2-'(l + h 2 )[l - P(h)\ - h(2*r m e- hi '\ 

4. When /i approaches - oo, K(r, ») approaches one and ./(r, /i) is 
approximately h + r. 

As t becomes very large pi(r, 7) approaches the value of N, which is 

(2tt)- 1 /3 1/2 exp (-/ 2 /2). (01) 

This follows from considerations of the type which led to (5). It may 
also be derived from (47), (48) and the fact thatm, m! , m" all approach 
zero as t becomes large. In this case M 22 approaches j8, and r and h 
become zero. 

When I is zero, equation (47) reduces to a result given in Ref. 2, 

namely 

Pl (r, 0) = AfdT w (l - m 2 r 3/2 [r(7r - cos" 1 r) + s ]/27r, (02) 

where < cos - r < tt. 

Expressions for the various parameters entering the expression for 
Pi(t, I) when t is small are given in Appendix II. Combining the results 
of Appendix II with (47), we obtain 

lh (r, I) W (r/J/4/3) exp (- lV/3/8)./(l, /i) 

/i^ -IPB~ U - 

where t is assumed to be small and /? is given by one of equations (125). 
There are two cases to be considered. In the first one, I remains fixed 
as t approaches zero. It then follows that pi(r, I) approaches zero 
linearly with t according to 

Pl (r,I) tt(rB/mJ(.l,h) (04) 

where .7(1, h) and h are given by (00) and (03). In the second case, I 
is positive and increases as t decreases so as to keep the order of tI at 
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unity. Then h approaches — a> and J(l, h) is approximately h 2 , with 
the result that 

Pl (r, /) » G3r7 2 /4) exp (-/V/3/8). (65) 

This expression may also be obtained from (47) by first letting 7 become 
large (assuming m! to be negative so that h becomes large and negative) 
and then letting r become small. 

When the right-hand side of (65) is integrated from t = to t = a> 
the result is unity. Since the same is true of the exact probability density 
p(r, I) [to which pi(r, 7) is a first approximation], and since pi(j, I) ^ 
p(r, I), we are led to expect that 

p(r, /) -* (/3r/ 2 /4) exp (-/V/3/8) (66) 

when I becomes large and positive (and rl remains fixed). A formal 
proof of (66) has been constructed on the assumptions (i) that p(r, I) 
ultimately decreases exponentially as t becomes large, and (ii) that m' 
remains negative for all positive values of r. The second assumption 
seems to be required by the method of proof rather than by nature. We 
shall not give the proof here because of its length. Instead, we shall 
outline a method by which (66) may be obtained directly.* This method 
is due to D. S. Palmer and grew out of correspondence with him. 

When I is very large, the intervals during which I{t) > / are very 
short. The length of an interval starting at t = U is, to a second degree 
Taylor approximation, 

I = 2I'(ti)/[-I''(t 1 )). 

The joint probability density of I(t) and its first two derivatives is 

where B is given by (125). Since I" is a random variable whose mean 
value — /37 is large compared to its standard deviation, we may put 
I = 2/'(/ 1 )//37. 

The chance that I(t) will pass upwards through the value I in the 
interval t\ , h -\- dt with a slope between /' and I' + dl' is equal to the 
chance that, at time tx , I{t) lies between I and I — I'dt and has a slope 
in /', /' + dl'. This chance is 



(I'dtKdl') T 7 2 7 

ex P - TT ~ 



'--[ 



27.V/3 L 2 2/3_ 



* Another proof has subsequently been devised by Mark Kac and D. Slepian. 
Also see the end of Appendix III. 
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Expressing /' in terms of I shows that the chance of I(t) passing up 
through 7 in t x , h + dt and starting an 1(1) > I interval whose length 
lies between I and I + dl is 

\j3 m fl dl dt/8ir] exp I" - ~ - ^-"1 . (67) 

Our goal of obtaining the probability density (65) for the lengths of 
the I(t) > 7 intervals, 7 very large, is now within reach. All that remains 
to be done is to divide (67) by the chance \p m dl/2ir] exp [-I 2 /2] that 
I(t) will pass upwards across 7 in ti , t\ + dt. 

In much of the following work it is convenient to use the variable u, 
defined by 

u = r/t, t = 7r/T 1/2 [l - P(/)] exp (/ 2 /2), (68) 

where the expression for the average duration t is obtained by setting 
b = 1 in (6) . The probability density f(u, I) of u is related to the proba- 
bility density p(r, I) of t by 

/(u, /) = tv(ul, I), (69) 

and the cumulative probability 

F(u,I) = J f(u,I)du (70) 

has already appeared, for the special w(f) given by (16), in Section III. 
The first approximation to /(it, /) corresponding to pi(r, /) is 

/i(«, /) = t Pl (ut, I). (71) 

When u becomes very large /i(u, 7) approaches the value [1 — P(7)]/2. 
The function /(m, I) has the following limiting forms: 

/(«, /) - (ut 2 B/^)J(l, h), u -» 

f(u, I) -♦ (ttu/2) exp (-mV/4), 7 -> co (72) 

/(m, 7) -> exp (-«), 7 -» - co, u > 0. 

The first of these follows from (64), the second from (66) and the second 
of equations (6) and the third expresses the exponential behavior which 
sets in after f(u, I) has passed through its maximum value (which 
occurs close to u — when 7 is large and negative). The exponential 
law is followed because most of the intervals are so long that the proba- 
bility of an interval of length r ending in t, t -\- dr is independent of r. 
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VI. APPROXIMATIONS FOR LONG INTERVALS 

It is generally assumed that, for large values of r, p(r, T) decreases 
exponentially. This assumption is supported by the following informal 
argument, which is much the same as the one used at the end of Section 
V for the case I = — °° . When T(t) remains greater than / for an inter- 
val long in comparison with the effective range of the autocorrelation, 
1(f) "forgets" when it entered the region I(t) > I. After this amnesia 
occurs it wanders about in the region in a fashion independent of its 
entrance time, until it happens to drop below the value / and thus end 
the interval. The exponential decrease is obtained when it is assumed 
that the chance of dropping below the value / in the interval t, t + At 
depends only on At and not on how long I{t) has been in the region. 

However, no one seems to have given an exact and usable asymptotic 
formula for this exponential decrease, although progress in this direction 
has been made by Kuznetsov, Stratonovich and Tikhonov. 9 When 
applied to our problem, their result gives F(u, I) as an exponential 
function of an infinite series whose nth term is a multiple integral of 
order n — 1. When t becomes large each term approaches the form 
a n T + b„ , where a n and b« are independent of r. Here we shall be con- 
tent to connect, by eye, the curve showing our estimates of p(r, /) for 
small values of t (obtained from pi(r, /)) with an exponential curve 
chosen so as to (i) make the area under p(r, /) unity and (ii) give the 
correct average value for r. It is hoped that this will give an estimate 
of some use in engineering applications, even though it may not con- 
tribute a great deal to our understanding of the behavior of p(r, I) for 
very large values of r. 

We shall work with u = t/1, its probability density f(u, I), denned 
by (69), and the approximation /i (u, I) denned by (71). The procedure 
consists of calculating 

F^uJ) = 1 - f fi(v,I)dv (73) 

Jo 

G x {u,I) = 1 - f" F 1 (v,I) dv (74) 

•><) 

by numerical integration and plotting F x versus (?i as shown by the solid 
line in Fig. 6. The curve starts at the point (1, 1), corresponding to 
u = 0, and curves in towards the origin as u increases, eventually cross- 
ing the G axis as shown. With the (Gi , F x ) curve as a guide, we draw in 
the (G, F) curve shown by the dashes. This curve is drawn so as to 
coincide with the solid curve in the region around (1, 1) and then depart 
from it in a "reasonable" manner so as to pass through the origin. Once 
the ((?, F) curve is determined, the relation between G and u may be 



W NATION OF FADES IN RADIO TRANSMISSION 



607 




Fig. 6 — Method used to estimate probability distribution l''(u, I) for large 
values of n . These curves are for 1=0 and the normal-law power spectrum (.16). 



determined by numerical integration from 

f 1 clG 

11 = L m 



(75) 



and from this the cumulative probability F(u, I) may be obtained as a 
function of a. 

The justification of this procedure follows from the fact that func- 
tions F(u, I), G(u, I), computed from /(u, I) (assumed known for the 
moment) by omitting the subscripts in (73) and (74), define a curve on 
Fig. (5 which starts at (1, 1) for u = and reaches the origin at u = oo. 
Thus 



F(oo,/) = 1-1 f(v,I)dv = 0, 

Jo 

G(x,/) = 1 - ( F(v,I)dv 

Jo 

= 1 - I vf(v,I)dv = 0, 
Jo 



(70) 
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( ti\ = ^ 

MT ' "' 47T(^ BB ) 1 /W 22 iV33Z) 3 

fl + a - b - c a - be f . . x ~| 

L r+~ a w*~ {ai + a2 + a:i _,r) J 

where AT,-,- is the cof actor of the element ij in the determinant 
1 
B 



m) 



N = 



n 



in —m 
in —m 
n —n 



m m n 

" i / 

—m —m —n 

B ri 

1 n 

n n 1 



with ii and n' being the values of m and m' computed for r/2. The other 
parameters are 

a - -N 23 (N 22 N 33 y m , b = -N 3b (N 33 N 66 y m , 

c = N 26 (N 22 N 66 r m , 

,_i a — be 



«i = cot 



0,1/2 ' 



2> 3 = 



1 


c 


b 


c 


1 


a 


b 


a 


1 



where < a Y < 7r and a 2 and a ;i are obtained from a x by a cyclic per- 
mutation of a, b, c. It turns out that b = c and a<> = a 3 . From the de- 
terminant AT it may be shown that 

ATa = N 33 = BA + B - m'\l - n 2 ), 
N, 3 - mM + B - m'V - m), 
A^25 = -N 3b = [n'(l + m) - nm'\C, 
N 56 = [03 + »*)(] + m) - m'']C, 
N = (N 22 2 - N 23 2 )/A, 
where 

A = (1 - m)(l + w - 2rf), B = n'(l - m)[2m'ii - n'(l + m)], 
C = (B - m"){\ - m) - m'\ 
The second item is concerned with the approximate integral equation 



DURATION OF FADES IN RADIO TRANSMISSION 611 

Pl ( Tj /) = p( T> I) + f p(»/, 7)«i(t - //, /) dy (84) 

where g x (T, /) dr is the chance that I(t) will pass downward through 
level 7 in (t, t + dr), given that /(/) passes downward through level 
I at r = 0. This is an extension of an equation given by D. S. Palmer. 
When 1 = 0, (84) becomes McFadden's integral equation mentioned 
in the discussion following equation (19). Equation (84) is based on the 
assumption that the length of an /(/) > I interval and the distance 
between the following downward crossings are independent. Changing 
the upper limit of integration for // from co to - » on the right-hand 
side of (38) gives an expression for q x {r, I) which leads to 



gi(r, /) = pi (r, /) + \m*#- 1, \\ - ™r 3/2 ex P \L - j-qr^]} 

J * exp(-A72)-^4^ [i-P(/01 



(85) 



in which the various parameters have the same meaning as in (47). 
When 1 = 0, the results (58) lead to 

(Ji (t, 0) = Ma0~ u \l - w 2 )~ 3/2 [s - r cos" 1 r]/2w. 

VII. DISTRIBUTION OF INTERVAL LENGTHS FOR ENVELOPE 

Here we study the envelope /?(/) of /(/) in much the same way as I{t) 
was studied in Section V. 1(1) is now taken to be a narrow-band Gaussian 
noise current whose power spectrum w(J) is assumed to be symmetrical 
about the midband frequency /„ . We shall take p{r, R) to be the proba- 
bility density of the interval lengths during which R(t) remains less 
than R. Strictly speaking, we should put some sort of a subscript, such 
as the /• used on F in Section III, on p(r, R) to distinguish it from the 
p(r, /) of Section V. However, we shall not do so because the chance of 
confusion is small. 

The first approximation to p(r, R) is pi(r, R), where pi(r, R.) dr is the 
chance that R(l) will pass upwards through level R in (/•> , U + dr) given 
that R{t) has passed downward through R at t Y , r being the difference 
k — h. Reasoning similar to that leading to (38) gives 

P>(r,R) = -feY R-V*'* [ dRl 

W °° k (86) 

f dRt' Ri R% p(R,Ri ,R* ,R) 

Ja 
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where now b = 1 and is the p r of (8). The function p(R t , Ri', R 2 ', R 2 ) 
is the probability density of R(h), R'(h), R'(t 2 ), R(t 2 ) and the primes 
denote time derivatives. 

Let I c (t) and I, if) denote the "in-phase" and "quadrature" compo- 
nents associated with R(t). Let 7 c i , lei, la', I a stand for the values of 
I c (t) and I J (J) at h and t 2 , and similarly for I e (t). Because of the sym- 
metry of w(f) about /o , the I c 's are independent of the I„'s and their 
joint distribution is 

p(Ici , In, U, IrtMLi , Li', h2', /. 2 ), (87) 

where these two functions may be obtained from (39) by adding the 
subscripts c and s, respectively, to the TVs and TVs. The /3 and m appear- 
ing in the expression (40) for M are now given by 

m = / w (/)cos 2ir(/ - / )t df, 

j ° (88) 

(* = [-m") T= , 

instead of (41). These statements follow from results given in Appendix 
II of Ref. 3. 

When we introduce the polar coordinates 



(89) 



(90) 



Id = R\ COS 0i lc2 = Ri COS 02 

I,i = Rx sin 0i 1,2 = R 2 sin 2 
and their associated time derivatives, and then set 

<p = 0i — 02 , c = cos <p, s = sin <p, 
i?i = R 2 = R, 
the joint distribution (87) goes into 

nr~ exp {" 2i7 [(2ilfu + cMu)r2 + 2(Mu + cMu ^ RRi ' - RR *') 

+ JtfnW" + ft'") + 2cM 2 ,R,'R 2 ' 

, 22 , (91) 

4- MnR'W 4- 2 ' ) + 2cM 23 R 2 d 1 '8 2 ' 

+ 2sM n R(R l 'd 2 ' - fl 2 '0i') - 2sM u R\d 2 ' + 0/)]}. 

The quantity piRfo'fa'ffldRidRdRidRz is obtained by multi- 
plying (91) by 

<tf c i • • • tf/ s2 = R A dRi dR 2 dRi' dR 2 d6 Y dd 2 dB x ' dd 2 ' 

and integrating 0i , 2 from to 2t, and 0/, 2 ' from — <» to + °o . One 
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of the integrations with respect to the 0's may be performed. We write 
the other as an integration with respect to <p. When the integrations 
with respect to Q\ and d 2 are performed p(R, R\ , R 2 , R) is found to be 



\2tt/ h 



in 

2 Jx,r 2\-l/2^ 



(M 2 i - c 2 Mn 2 r'V d<p (92) 



where A contains terms of the first and second degree in Ri, R 2 '. Inci- 
dentally, the linear terms enter only in the combination Ri — R 2 . 
Upon setting Ri, R 2 equal to Xi — a\ , x 2 — a 2 , respectively, and 
choosing a x , a 2 so as to make the coefficients of Xi , x 2 in the resulting 
expression for A vanish we find 

a 2 = —ai , (93) 

CMW - M 23 2 ) ai = R[MnMn + M»Mu + c[M 23 M n + M 22 M 2i ]}. 

The expressions given in (123) show that both sides of the equation for 
ai contain the factor M, which may be divided out with the result 

= Rm'jc - m) m 

1 — m- 

This choice of ax and a 2 and the results of Appendix I enable us to reduce 
(after a considerable amount of algebra) the expression for the exponent 
A in (92) to 

(1 - m) r , f ( t , ax , „„,, „ „ 1 /g 2 (l - mc) , Qr , 
M/2SOT1 + x 2 ) + ZcAl 23 x h v 2 \ — — : — . (95; 



2(A/ 22 2 - «) 1 -w . --, . 1 - m 2 

When (92) is put in the expression (86) for pi(r, R) the result is a 
triple integral. The expression (95) suggests that we replace the variables 
Ri, R2 by x, y defined by [see (46)] 

x = -.r,[(l - m 2 )/M 22 } v \ y = .r 2 [(l - m*)/Mn] m , 

(96) 
p = CM23/M22 ■ 

Then (86) goes into 

t m RM 22 e R *' 2 f 2 * , n T ft 2 (l - wic) 1 . , Q7 , 

Pl(r ' B) = (2*py<Ki - my J. J(p ' fc) exp L" 1 - «■ J * (97) 

where «/(p, /c) is obtained from (48) by putting p, k for r, h and 
k = -ai[(l - m 2 )/M n ) m 

Rm'(m -c)(\- m 2 Y /2 (98) 



1 — m- 



( I - m 2 V /2 
\ M 22 ) ' 
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in which c stands for cos (p. Here M& and M& are obtained by setting 
the values (88) for m and /3 in the appropriate expressions given in 
(120). 

When t becomes very large, k and p approach zero and we see that 
Pi(«, R) = R(p/2r) m exp (-R 2 /2) as it should according to (9). 
When t approaches zero, with R fixed, the results of Appendix II enable 
us to show that 

Pl (r, R) -> -^L f J[l, RfiB- m (l - x^e-*'*" 2 dx 

2/3v2tt J o 

where ./(I, h) is given by ((50). 

The integral (97) must be evaluated numerically. However, a special 
case of interest in fading is obtained when R approaches zero. In this 
case, Pi(t, R) tends to depend on r and R only through the ratio t/R. 
When we let ^ and R both approach zero in such a way as to keep t/R 
fixed, the results of Appendix II show that 

il/ 22 B R\l - me) _ R\l - c) 



(1 - m 2 ) 2 4/3 ' 1 - m 2 0r 2 ' (99) 

A; — > —22(1 — c)2B~ m T~\ 

Therefore, k approaches — 00 except when <p is zero, and J(p, k) is 
approximately l>~ except in a negligibly small region around ip = 0. It 
is convenient to introduce the variable u = t/1, which, in the present 
case, becomes, with the help of the approximation (11) for I, 

—1(0"- (I00) 

The function of u corresponding to pi(r, R) is 

/,(», R) = 1 Pi (t, R), (101) 

and it follows from (99), (100) and (101) that as R approaches zero, 

Mu, R)-> Z - I'' (1 - c) 2 exp [- 2(1 - c)] d<p, 

2, Jo (102) 

z = 2/(ttm 2 ). 

When we write (1 — c) 2 as i — 2c(l — c) and integrate the s 2 por- 
tion by parts, we find that (102) may be expressed as a derivative: 

h(u, R) -> - i- — I c exp [- 27r'<r 2 (l - c)] d<p 

(III TTU Jo 

(103) 

d 



du 



u J 
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where I\{z) is a Bessel function of imaginary argument. It may be shown 
that the expression under the differentiation sign is unity when u = 
and is zero when u = cc . Hence, an argument similar to the one follow- 
ing equation (05) suggests that as R approaches zero the approximation 
Pi(r, R) approaches the actual probability density p(t, R). Then/i(w, R) 
approaches /(w, R) and (103) is the limiting form of both. We also have 



F(u, R)->- /i(2/tm 2 ) exp (- 2 /mi 2 ) 
u 



(104) 



u - 


->o, 


u - 


-* °°, 


f(u, R) - 


-> 3iru/8, 


/(«, R) - 


-* G/wu , 


F(u, R) - 


1 1 2 

-» 1 - fsiru , 


F(u, R) - 


-> 2/W. 



where F(m, 7?) is the probability that an interval during which R(l) < R 
has a (normalized) length greater than u. F(u, R) is the function de- 
noted by F r (u, R) in Section III where equation (104) appears as equa- 
tion (23). When u is very small or very large we have 



(105) 



The indefinite integral of the right-hand side of (104) is a constant plus 
<T-[/o(z) + hi?)], * = 2/W. (100) 

Equation (97) refers to the intervals for which R(t) < R. If we are 
interested in the intervals for which R(t) > R, the analogue p*(t, R) 
of 7>i(r, R) is obtained from (80) by changing the signs of the »'s in 
the limits of integration. Instead of J(p, k) in the integrand of (97) we 
get a double integral similar to (48) except that the upper limits 
of integration are — oo's instead of w's. This changes the J(p, k) 
to J(p, —k). Upon using (50), we obtain 

Pi*(t,R) = lh (r,R) 
««/2 



-A.- 2 /2 



BMvf ' 2 r\ ( , j^ p(] , 2k 

(2irj8) 1/2 (l - m-)- Jo L V^r 

R\l - me) 



(107) 



exp 



1 — m 2 



dtp 



which is the analogue of (57). When R becomes very large and r small 
in such a way to keep Rt fixed, pi(r, R) approaches zero from physical 
considerations. We also have near <p — 



#; _ fl 2 (l - mc) 
2 1 - m 2 



0ft 2 



2/St 2 ' 



(108) 



BjSB -1 ' 2 - R<p 2 T- 2 B- m . 
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The first line in (108) shows that most of the contribution to the inte- 
gral in (107) comes from the region where <p is 0(t/R). In this region, 
the second term in the approximation for k is negligible in comparison 
with the first, so k is large and positive and hence P(k) is nearly unity. 
Thus when R is large (107) becomes 



Pi*(t,R) -> — — e 



(109) 



which is of the same form as (66). Therefore, as we might expect, the 
interval lengths for large values of R are distributed in much the same 
way as are the interval lengths of I(t) for large values of /. 



VIII. COMPUTATIONS 



The computations for the curves shown in Figs. 1 and 2 were based 
on either the power spectrum (16) or (20), depending on whether I(t) 
or R(t) was being considered. In both cases the bandwidth parameter a 
was taken to be l/(27r), so that both /3 and /3 r were unity and 



m = exp (— t /2). 



(110) 



This autocorrelation function was also used by Palmer in his computa- 
tions. 

Fig. 7 shows curves of the first approximation /i(w, /) to the probabil- 
ity density f(u, I) for the lengths of the intervals during which I{t) 
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Fig. 7 — Graphs of the first approximation fi(u, I) when the power spectrum 
has the normal-law form (16). The rms value of l(t) is one. 
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2.2 2.4 



Fig. 8 — Graphs of the first approximation F i(u, I) to the interval length 
distribution F(u, I). Estimates of /*'(?/, /) are shown in Fig. 1. 

exceeds I. These curves were calculated from (47), (71) and tables 8 
giving values of K(r, h). Fig. 8 gives curves of Fi(u, I) computed from 
/i(w, 7) by means of (73). Fi(u, I) is the first approximation to the cumu- 
lative distribution F(u, I), estimates of which are shown hi Fig. 1. One 
way of estimating F(u, I) from F x (u, I) is to connect the F^u, I) curve 
smoothly to a curve which approaches a straight line on the semi-log 
coordinates of Fig. 8, the slope and position of the straight line being 
chosen so as to make 



[ F{u, I) du = 1. 
Jo 



(HI) 



This relation follows from the fact that the average value of u is unity. 
The linear behavior ensures that F(u, I) decreases exponentially as u 
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Fig. 9 — Graphs of (Gi, I'\) curves and estimates of the corresponding true 
(G, F) curves for 1(1). 



becomes large. However, here F(u, I) is estimated from Fi(u, I) by 
the method outlined in Section VI. Fig. 9 shows (G, F) curves of the 
type illustrated in Fig. G. The estimates of F(u, I) given by Fig. 1 were 
obtained by connecting the curves of Fig. 9 to lines passing through 
the origin in a "smooth and reasonable" manner. The relation between 
u and F(u, I) was determined by numerical integration from (75). 

As a guide to what "smooth and reasonable" means, the case 1 = 
was studied more completely with the help of the second approximation 
(83) and the computations made by Palmer 5 from his approximate 
integral equation. Some of this information is shown in Fig. 10. Palmer's 
results lead to a (G, F) curve (not shown) which lies slightly above the 
second approximation curve when G lies between 0.25 and 0.55. As G 
decreases below 0.25 the separation between the two curves increases 
and the curve from Palmer's results coincides with our estimate of the 
true (G, F) curve. 

Fig. 11 shows curves of the first approximation f\{u, R) to the proba- 
bility density /(«, R) for the lengths of the intervals during which the 
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Fig. 10 — First and second approximation (G, F) curves for / = 0. 



1,2 



0.2 -i 



R=2 


i.ieJOcs' 


[0.5 

P^o 
















7 / /A 














0.865 




























\ 




















\\ 








0.5 


/ // // 


^\\ 


\ 




0.393 




/// 




































_ 



0.2 0.4 0.6 



1.0 

u 



1.2 



1.6 1.8 2.0 



Fig. 11 — Graphs of the first approximation fi(u, R) when the power spectrum 
has the normal-law form (20). The median value of Nil) is 1.18. 
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Fig. 12 — Graphs of the (G u F t ) curves for R(l). 



envelope R(t) is less than R. These curves are the analogues of the 
curves for 1(f) shown in Fig. 7. Likewise, the (G, F) curves shown in 
Fig. 12 are the analogues of the curves of Fig. 9, and the estimates of 
F r (u, R) given in Table III were obtained from them. The functions 
fi(u, R), defined by (101), were computed by numerical integration from 
(97). 



IX. COMMENTS ON A MARKOFF PROCESS 

The foregoing work supposes that the autocorrelation function m of 
1(1) is such that its second derivative with respect to r exists at r = 0. 
This is not true when w(f) has the form 4/[l + (2ir/) ]. In this case m 
is exp (— | t \) and m" is infinite at r = 0. I(t) is now a "Markoff proc- 
ess" and a formal calculation of the average distance between successive 
zeros gives l = 0. 

Although the distribution function for the intervals between zeros 
does not exist in the usual sense, some related information is available. 
Most of it is due to A. J. F. Siegert. 12 Here we shall list several of his 
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results and comment on their relation to the work of this paper.* For 
the sake of simplicity we take /"(/) = 1 and m = exp (— | t |). 

1. The probability W(I, r) dr that 7(0 reaches the value for the 
first time between r and r + dr ("first passage time") starting from 7 
at time t = is 

W(I, r) dr = (2/V) I/2 / exp (-fz 2 /2)(-dz), 

z = e -(i - e - 2 T m , -dz = eV dr. (112) 

Wang and Uhlenbeck 13 obtained this result by solving an associated 
diffusion problem (by a method first used by Smoluchowski) in which 
the I = axis corresponds to an absorbing barrier. Another derivation 
of (112) has been given by Siegert, 14 who obtains it as a solution of an 
integral equation. 

2. The chance that 7(0 will go through zero for the first time in t, 
t + dr given that 7(0) > is 

Q(r) dr = (2A)<T T (1 - e~ 2 T m dr. (113) 

By symmetry, the condition 1(0) > may be replaced by 1(0) ^ 0. 
Q(r) dr is also the probability that the interval between an instant se- 
lected at random and the next following zero of 7(0 has a length / lying 
between t and r + dr. 

3. The chance that I exceeds r seconds is 



f 



Q(t)dt = -sin-'e-' (114) 



where the arc lies between and ir/2. 

4. Let 7(0 be observed over a long period of time T. The expected 
number of spacings between successive zeros whose lengths exceed r 
seconds is 

N t (t) = Q(r)T. (115) 

Furthermore, the probability that an instant selected at random will 
fall in a space whose length exceeds r is 

tQ(t) + (2/tt) sin -1 e~ T . (11G) 

Items 2, 3, and 4 are due to Siegert 12 and were obtained by analysis 
similar to that given in his paper. 14 David Slepian has pointed out, in 



* I am indebted to Mr. Siegert for a copy of his report, which contains the 
results listed here. 
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some unpublished work, that (114) may also be obtained from the 
probability 

P++ = i + 1 sin" 1 <T (117) 

4 zt 

that both 7(0) and 7(t) be positive. His derivation runs as follows. The 
probability (114) is twice the probability that I(t) be positive over an 
interval of length r selected at random. This last probability is P ++ — 
P+o+ where P+o+ is the chance that 7(0) > 0, 7(r) > and 1(1) is zero 
at one or more points in < t < r. Since the process is Markoffian and 
symmetrical about zero, consideration of the behavior of I(t) between 
its first zero (if it has one) and r shows that P+o+ is equal to P+o_ , 
i.e. to P + _ . The sum of P ++ and P + _ is the probability that 7(0) be 
positive, which is \. We then have 

2(7^ ++ - P +0+ ) = 2(P ++ - P + _) = 47 J ++ - 1 = - sin' 1 <T T 

T 

and this gives (114). The expression (117) for P ++ is readily obtained 
by integrating the joint probability density. 

5. It should be possible to obtain (114) by letting n approach infinity 
in the expression for two times the probability that 7(0), 7(r/n), • • •, 
I{nr/n) all be positive. The expression which should approach (114) is 

a,-cfi)/i (1 _ a y/ 2 r dXi ... r dXn+1 . 

Jo Jo 

exp [- .rf - -W - (1 + a)(xi + • • • + x,;) ( 118 ) 

4- 2a(xiX 2 4- x. 2 x 3 4- • • • 4- x n x„ + i)] 

where a denotes exp (— r/n). 

This method of obtaining (1 14) is appealing because of its directness. 
Unfortunately, it seems that no one has yet succeeded in carrying out 
the limiting process. Some idea of how the limit is approached may be 
obtained from Table IV. 

The values shown in the columns for n = 1 and n = 2 were computed 
from the corresponding known expressions for (118), namely 

i 1 — 1 -t or> 

1 cos e = 2r ++ , 

IT 

1 - (£)(2 cos -1 e" T/2 4- cos -1 e~ T ). 
When t is small these behave like 1 — 0.45 \/t and 1 — 0.54 \/r 
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Table IV — Values of Expressions (118) and (114) 







Value of (118) 






Value of (114) 


T 


e~ T 


» = 1 


M = 2 


n = 3 


M = 4 





1 


1 


1 


1 


1 


1 


0.05 


0.951 


0.900 


0.879 


0.868 


0.861 


0.800 


0.10 


0.905 


0.860 


0.830 


0.814 


0.804 


0.720 


0.15 


0.861 


0.830 


0.793 


0.774 


0.761 


0.660 


0.20 


0.819 


0.805 


0.763 


0.742 


0.726 


0.611 


0.40 


0.670 


0.734 


0.672 


0.641 


0.620 


0.474 


oo 





0.500 


0.250 


0.125 


0.0625 






respectively, while (114) becomes 1 - 0.90 Vr- When n = 4 the five- 
fold integral may be reduced to 



2x- 1/2 (l - a 2 ) 172 f dx 3 [K(r, -h, -k)} 2 exp [-(1 - a*)xf\ 

Ja 



(119) 



where K(r, h, k) is the d/N tabulated in Ref. 8, and 



2\-l/2 



.1/2 



k = rh. 



r = a(l 4- aT l '\ h = 2 u -ax 3 , 

When a is near unity it is helpful to set 

K(r, -h, -k) = 1 - L, 

L=l- K{r t K *) - hP(h) - $P&). 

The contribution of the 1 in 1 - 2L 4- Is to (119) is unity and the 
contribution of -2L + If may be obtained by numerical integration. 
Replacing the square of K{r, -h, -k) in (119) by 

K(r, -h, -k)[l + P(h)]/2 

gives an expression for (118) when n = 3. 

6. The probabilities stated in Items 2, 3, and 4 have already appeared 
in Section VI for the well-behaved cases in which t exists. Thus, from 
(79) the analogue of Q(t) dr is P(t, I) dr/t, or F(u, I) du, and the analogue 
of (114) is G(u, I). The probability (81) that a point chosen at random 
falls in an interval longer than t goes into the corresponding expression 
(116). 

These analogies lead us to regard Q(t) and its integral (114) as being 
irregular versions of F and G (with t playing the role of u) for which 
a (G, F) curve may be plotted with the help of expressions (113) and 
(114). When this is done, it is seen that F approaches infinity as G 
approaches unity. However, the other parts of the curve behave in 
much the same manner as in the earlier cases. 
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APPENDIX I 

Minors of the Determinant M 

Here we list the values of the minors of the determinant M denned 
by equation (40) and state some of the relations between them. 

In = M u = /3 2 - to" 2 - /3m' 2 , 

— Ma = M 3i = m'm" + mm' (3, 

Mis = -M 2i = m' 3 - m' /3 - mm'm", (120) 

M 22 = M 33 = /3(1 - to 2 ) - m'\ 

Mu = mm"* — m/3 2 — m' 2 m", 

M 23 = m"{\ - w 2 ) + toto' 2 , 

M = [(/3 + to")(1 + to) - to /2 ][(/3 - to")(1 - to) - m'\ (121) 

Mn + M u = M n + M 2i = M n - M 23 
P + to" to' 1 + to 

= (/3 - to")(1 - to) - to' 

M 22 4- M n = (1 - m)[(|8 + to")(1 4- to) - to' 2 ] 

M = (M n 4- M u )(l + to) - m'(M a 4- M u ). 

We also have relations of the form 

Mo 



2 - M^ = 


M22 M n 
M 32 M 33 


/\2+2+3+3 


1 TO 
TO 1 


M 






= (1 - w 2 )M, 


|- M n M 2i = 


M n Mi 3 
M 32 M 33 


— /•_•) 1+2+3+3 


-to' 

TO 1 


M 








= — to' 


mA 


f, 



(122) 



(123) 
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MvMvi + M22M24 = m'M , 

Mnitf M - Mj = OS - m /2 )M, 
M U M 2Z + il/ 13 2 = (ro /2 - mm")M, 
MuMsm + M 14 il/ 22 = (m" - /3m) M. 

APPENDIX II 

Expressions for Various Parameters When r Becomes Small 
In dealing with 1(1) we shall put 

b 2n = (27r) 2n f «>(/)/'" #, » = 0, 1, 2, • • • (124) 

and in particular 

6o =1, b 2 = 0, b 4 = 7, 

5 = 6 4 - W = 7 - P\ (125) 

C = b 2 b 6 - b 4 2 = 0b 6 - 7 2 . 

The results for I{t) may be carried over directly into those for the en- 
velope R(t) of a narrow-band noise current whose power spectrum is 
symmetrical about the midband frequency f . This is accomplished by 
replacing the definition (124) for b 2 „ by 

6 2 n = (2t) 2b [ w(f)(f - / ) 2n df. (126) 

For the normal-law power spectrum used in our computations 7 is 
equal to 3/3 2 and B becomes 2/3 2 . 

It follows from (41) that when t is small and the b's exist, 

/3t 2 7t 4 b 6 r 

m = 1 "2T + 4T""6T + ""' 

* , (127) 

m = —@t + • • • , 

OT " = = -p + . . . . 
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These expansions may be used to show that 

1 — m 



2 ' 



1 + m _>2-^- 



9, 



1 2 n 2 

1 — m —> Br — 



08 + m"){\ + m) - m' -* r 2 B, 
(8 - m")(l - m) - m'* -* r 6 C/144, 

il/ 23 -» i¥ 22 -> r 4 /?£/4, 
il/22 - il/23 -» r 6 C/72, 

M -> r 8 J5C/144, 
il# 22 2 - M 23 2 -» r ,0 /3BC/144, 
(1 - m 2 )/M 22 -^ 4/^r 2 . 



(128) 



APPENDIX III 



An Approximation to the Cumulative Probability P(t, I) 

In Section V we worked with a first approximation pi(r, /) to the 
probability density p(r, /) for the lengths of the I(t) > I intervals. 
Here we note a rather simple approximation Pi(t, I) to the probability 
P(t, I) that an interval selected at random from the universe of /(/) > / 
intervals will be longer than r. We take Pi(t, I) to be the limit, as dt 
approaches zero, of a fraction having the numerator 

[Probability that /(r) > / given that I(t) passes upward 

through the value / during the interval (0, dt)] 

and the denominator 

[Probability I(t) passes upward through / in (0, dt)]. 

By proceeding as in Section V it may be shown that 

Pi(t, D = I ~ \ P (*) + \ [l + P (y)l exp (-/,- 2 /2) (129) 
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where P(x) is the error integral (3), and now 

r = - m '(r 1,2 (l - m 2 )' 1 ' 2 , s 2 = 1 - r\ 

k = /[(l - m)/(l + m)] v \ 

with w and /3 the same as in Section V. The approximation given by 
Pxiu't, I) to F(u, I) is not as good as the Fi(u, I) given by (73) (which 
is based on the pi(r, /) of Section V), but we do have P\(ut, I) ^ 
F(u, I) ^ I'\{u, I). This may lead to another proof of (66). 

APPENDIX IV 

Fade Length Distributions Having a Large Percentage of Long Fades 

Some, but not all, of the small number of available experimental 
observations of fade lengths show a greater percentage of long fades 
than those shown in Figs. 2 and 3. Since our calculations are based on a 
normal-law power spectrum, it is of interest to know whether a different 
power spectrum may give a larger percentage of long fades. An example 
of such a power spectrum is studied in the first part of this appendix. 
The second part is concerned with a guess as to how the very long fades 
are distributed. 

The power spectrum 

.: ' *» = itw)" (130) 

which has been discussed in Section IX, is quite different from the nor- 
mal-law type. The power carried by the high frequencies in this spectrum 
is so huge that I = 0, and our method of investigating interval lengths 
breaks down. Our formulas will still apply, though, if the high frequen- 
cies are attenuated by some factor such as a A /[a + (2n/) ] . An attenua- 
tion of this type might either exist in nature or be produced by the 
measuring apparatus. Here we shall take 

m 8a\a + l)7(2a + 1) mn 

wU) = [i + wm + WW ( ° 

where, the numerator has been chosen so as to make the integral of w(f) 
from to ao equal to unity. 

If it were not for the considerable amount of additional work required, 
we would investigate the interval length distribution associated with 
the envelope of an I(t) having the power spectrum obtained from (131) 
by setting (/ - /o) in place of / and changing the 8 in the numerator to 
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u = r/t 

Fig. 13 — The probability F(u, I) that 1(1) > I for an interval lasting longer 
than t. 1(1) has the power spectrum given by equation (131) and t is ir(r l \/2a -+- 1. 

4. Instead, we shall apply the results of Section V to obtain the interval 
length distribution associated with an I(t) having the power spectrum 
(131), and rely on the similarity between the curves of Figs. 1 and 2. 

The general theory of Section VII shows that, as R approaches 0, the 
fade length distribution F r (w, R) of the envelope approaches the form 
F r (u, R -» 0) shown in Fig. 4. However, when a is very large it appears 
that this limiting form will have to be approached by going through 
shapes similar to those shown in Fig. 13. This indicates that the values of 
R will have to be very small indeed before F T (u, R) begins to look like the 
F r (u, R -+ 0) of Fig. 4. 

The autocorrelation function m, defined by (41), and its derivatives 
with respect to t are, assuming r > 0, 

= x — y[3a — 1 + ar(a — 1)], 

= -x + ya[2a + aria 2 - 1)], 

= x — ya[a + 1 + ar(a — 1)], 



m 



m' 



(132) 



>/;. 



where 



x = 2Aa 3 e 



y = A e 



A = l/(2a + l)(a - l) 2 . 
The expansion of m in powers of t up to the fifth is 

3 I r I 5 2(a + 1) J 



o 2 4 



5! 



2a 4- 1 



(133) 



(134) 
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in which 

(3 = a 2 /(2a +1), 7 = a (a + 2)/ (2a + 1), 

(135) 
B = y - 0- - 2a 3 (a 4- l) 2 /(2a + 1) . 

From /3 and 7 it follows that the average interval between the zeros of 
I(t) is ra~ l \/2a -f- 1, and the expected number of maxima per second 
is (2ir)~Va(o + 2). 

These expressions were used to compute pi(r, I) from (47), and the 
method outlined in Section VI was applied to obtain the curves shown 
in Fig. 13. These curves are for a = 10 and a = 100. Difficulty was en- 
countered, in the a = 100 case, in extrapolating the distributions for 
small values of t to obtain F(u, I) for the larger values of u. Conse- 
quently, these curves are not too reliable, but there seems to be enough 
truth in them to show that the percentage of long fades is larger than 
for the distributions of Fig. 1 pertaining to the normal-law power spec- 
trum. 

Irrespective of whether power spectra of the form (131) have any 
practical significance, it is interesting to speculate on how the corre- 
sponding interval length distribution F(u, I) behaves for large values of 
u. When a becomes large, I{t) has many more maxima than zeros (the 
average number of maxima between successive zeros approaches \/a/2) , 
and we are led to picture the power spectrum (131) as consisting of two 
parts. 

One part is the "low-frequency" portion which behaves like (130) and 
extends from zero up to a frequency at which the cutoff factor begins to 
operate. The other, the high-frequency portion, extends from this fre- 
quency to infinity. Most of the power is in the low-frequency part and 
we may regard it as producing most of the deviation of I(t) from zero. 
The high-frequency part produces small rapid fluctuations which are 
superposed on the more slowly changing low-frequency portion of I(t). 
Thus the high-frequency part produces the large number of maxima and 
the low-frequency part controls the drift towards or away from the value 
I(t) = 0. According to this picture some relation might be expected 
between F(n, 0) and Q(t) [given by (113)] when u is large. In fact, com- 
parison of (79) and (113) suggests that F(u, 0) du approaches Q(t) dr 
for the larger values of t: 

F(u, 0) -» l(2/v)s~ r (l - e~ 2 T m , r = id. (136) 

The curves shown in Fig. 13 indicate that (136) may not hold when 
a = 10, but may do so when a = 100. Evidently, if (136) is to be valid, it 
is necessary that the number of maxima between zeros be large. This is 



630 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1958 

also suggested by the original derivation of (112), which is based upon 
an analogy with the Brownian motion of a particle subject to an elastic 
force. 

If the relation (136) between F(u, 0) and Q(t) should turn out to be 
true, it would be helpful to have expressions for the generalizations 
Q(t, I) and Q r (r, R) of Q(t) corresponding to the "first passage times" 
down through the level /(/) = I and up through the level R(t) = R, 
respectively. These generalizations may be obtained by the method used 
by Wang and Uhlenbeck mentioned in Item 1 of Section IX. The deriva- 
tions are sketched in the following paragraphs. 

Let I(t) have the power spectrum (130) and let W(I, h , r)dr be the 
probability that I(t) passes down through the level I for the first time in 
the interval t, t + dr when it starts from the level Ii , which is greater 
than I, at time t = 0. Then W(I, 7, , t) dr is equal to (dP/dx) dr eval- 
uated at x = I and / = t, where P is the solution of the diffusion equa- 
tion 

*Z = il + x a J- + P (137) 

dt dx- dx 

which is zero at x = / for all time and is concentrated around x = h 
at time [i.e., equals 8(x - h) at time 1 = 0]. The probability Q(t, I) dr 
that if 1(0) > I then 7(0 passes down through the level I for the first 
time in r, t + dr is obtained from 

Q(r, I) = /" 6XP V2^" /2) W(I> II ' T) rf/l /P r ° b[/(0) > 7] - (138) 
A typical solution of (137) is 



- xVi D n (x) e'" 1 (139) 



where D u (x) is the parabolic cylinder function. The boundary condition 
at x = I gives the equation 

D n (I) = (140) 

to determine the eigenvalues n, , i — 1, 2, • ■ • . A formal application of 
the Sturm-Liouville theory shows that P may be expressed as a sum of 
terms of the type (139), with suitable coefficients, when n is replaced by 
Hi and the summation taken over i. If the typical eigenfunction is 
denoted by a,iD ni (x), the relation to determine the normalizing constant 
a,: is obtained by differentiating 

(n - m) f D„(x)D ni (x) dx = -DM)D ni '(I) (141) 
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with respect to n and then setting n = ??.,-. In (141), which may be 
obtained from the differential equation for D„(x), we have denoted 
dD n (x)/dx by D n '(x). 

It is found that 

W (7, 7i , t) = exp [(/f - 7 2 )/4] T.arD ni {h)D ni '{I) e~" iT , (142) 

1=1 

/o \-l/2 -' 2 I2 co t^ l,,s -»{T 

Q(t ' 7) = 2-U - P(/)] S n,-[-dZ) 1 ,(/)/dn] l .-. i ' (143) 

where P(7) is the error integral (3) and 

ar = D ni '{I) [-dD n (I)/dnU ni . (144) 

The integration (138) used to obtain (143) from (142) makes use of a 
result obtained by integrating a form of the differential equation for 

Dn(x). 

Siegert's result (115) now reads as follows: When I(t) is observed 
over a long period of time T, the expected number of I(t) > 7 intervals 
whose lengths exceed t seconds is 

NAt, I) = I 11 - P(I)]Q(t,7). (145) 

This may be seen by writing (79) as 

N r (r, I) (It , , 

Tprob [/(*)> / I = Q(T ' 7) dT ' 

When 1 = 0, the values of the n, are 1, 3, 5, • • • and it may be shown 
that (143) reduces to (113) as it should. As I becomes large and positive 
the n,'s increase in value, so that ni tt 7 2 /4 and the first few succeeding 
n/s differ from ii\ by terms of order I 213 . Furthermore, setting n = n^ 
and approximating D„_i(7) by (—1) dD„(I)/dn in the relation 

D„'(I) -r-^DM) -n7)„_ 1 (7) = 

suggests that when I is large and i not too large 

D n ' i {I)/n i [-dD ll {l)/dn], li ^\, 

and hence we may expect Q(t, I) to be of the form of 7 exp ( — 7 2 r/4) 
times a more slowly varying function of r and 7. 

In much the same way, expressions may be derived for W r (R, Ri , t) 
and Q t (t, R) corresponding to the envelope R(() of a Gaussian noise 
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current I(t) whose power spectrum is obtained from (130) by substitut- 
ing / — /o for / and replacing the 4 in the numerator by 2. Instead of 
(137) we must use 



dP 

dt 



» + (r-l)* + (l + l)p (146) 

dr 2 \ r ) dr \ r 2 / 



where P = 2irr P x P y now is the probability density or concentr ation 
corresponding to a narrow circular ring of radius r = \/x 2 + y 2 and 
P x , P y are, respectively, solutions of (137) and the corresponding equa- 
tion in y. 

The probability W r (R, Ri , r) dr that the envelope R(t) passes up 
through the value R for the first time in t, r + dr when it starts from 
Ri , less than R, at time t = is — (dP/dr) dr evaluated at r = R and 
t = t. Here P is the solution of (146) which is zero at the boundary 
r = R for all t and is 8(r — Ri) at time t = 0. 

The typical solution of (146) is 

■ 1/2 o- r2/ Va, r) e~ at (147) 



r e 



where a is a parameter similar to n in the earlier case, and 

ifa,r)-f^.^V.(-|;i;0 (148) 

becomes the typical orthogonal function for the range ^ r ^ R when 
a is one of the roots a t - , i = 1, 2, 3, ■ • ■ of the equation 

lFl (_| ;1; |) = . (149) 

For very large values of R, a t approaches 0, 2, 4, • • • and the confluent 
hypergeometric function in (148) becomes the Laguerre polynomial 
Li_i(r 2 /2). For very small values of R, v(a { , r) approaches 

Vr Jo(r\/ai + 1) 

and Rs/oci + 1 becomes the ith root of the Bessel function J (x). Some 
rough work indicates that «i is approximately 5 for R = 1 and 3.25 for 
R = 1.18. When R is very small «i +1 is approximately 5.78/R 2 . 
It is found that 

P = (r/Bd 11 * exp [(itf - r 2 )/4] E a«V«< , fl,M«< , r) C* 1 (150) 



1=1 
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where the normalizing constant a,- is given by 

= 1. 



2 V dvia, ft) "| r dv(ai,r) ~\ 
°'" _ da Ja= a , L dr J 



= R 



When the expression for W,(R, fti , r) obtained by differentiating P is 
set in the equation corresponding to (138), we obtain 

<Kr.g)- * T ±— [ £Hm7n (151) 

exp [Br/2] — 1 .=1 a,- ldv(a, ft)/daja=ai 

where Q r (r, R) dr is the probability that R(t) will pass up through the 
level ft for the first time in r, r + dr, given that ft(0) < ft. 
In the limit as ft becomes small, (151) reduces to 



.—air 



Q t (t, ft) -> 4ft 2 g e , (152) 

ft\/«i — * frth root of Jo(^)]- 

When (152) is integrated with respect to r between the limits and °o 
the result is 1, as it should be. The result (152) has probably been given 
before, since it is closely related to Brownian motion of a particle and to 
the first passage time across a circle, but a suitable reference has not 
been found by the author. 

When ft is small so that prob [R(t) < ft] -> R 2 /2, the relation (152) 
and the analogue of (145) for R(t) lead to 

00 

N t (t, ft) -> 2T E C" - " (153) 

for the expected number, in the long time T, of R(t) < ft intervals whose 
lengths exceed t seconds. 

All of the discussion between equations (136) and (153) pertains to 
the power spectrum (130) or its analogue for R(t). In the comparison of 
theory and experiment given in Section IV, it is convenient to deal 
instead with 

•W-FTOJi- (154) 

The analogue of (154) for R(t) has already appeared as equation (24). 
Results for (154) pertaining to duration of intervals, time averages and 
so on may be obtained from the corresponding results for (130) by re- 
placing t by kt. Thus, for example, the probability Q r (r, ft) dr goes into 
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kQrikr, R) cIt and the number of R(t) < R intervals longer than t to be 
expected in time T, when R is small, is 

=0 

N t (t, R) -> 2kT £ e~ aikT , (155) 

where a,- is still given by the second relation in (152). Siegert's expression 
(115) for the expected number of spacings between successive zeros 
whose lengths exceed r seconds is now 

N t (t) = (2kT/w) e~ kT (l - e- lk T u \ 

Incidentally, this may possibly lead to an estimate of A- when one is 
trying to measure a power spectrum which insists on increasing as 
\/f at the lowest frequencies which can be conveniently measured. The 
problem of estimating A; from N T {r) has been discussed by Siegert. 12 

If the "low-frequency" portion of a power spectrum behaves like 
(154) instead of (130) the remarks in the preceding paragraph show that 
the analogue of (136) is 

F(u, 0) -+ kt(2/ir) e~ kT {\ - e ~ u ' T )~ m (156) 

where u = r/i is assumed to be very large. The analogous surmise for 
F r iu, R) is that 

F r (u, R) -> klQ r {kr, R) (157) 

where the low-frequency behavior of the power spectrum is given by 
(24). When R is small, this and (152) suggest that 

F r (u, R) -> 4klR~ 2 £ e~ aikT (158) 

i=i 

for the probability that a long fade below the level R will last longer than 
t seconds. Multiplying both sides of (158) by N R and using 

Nnt = prob [R(t) < R] -* R?/2 

shows that the expected number of fades below level R, in unit time, 
whose lengths exceed r is 

N R F T (u, R) -> 2k E e- aikT (159) 

i=i 

where a,- is related to the zeros of J (x) as indicated by (152). The rela- 
tion (159) is to be expected in view of (155). 



DURATION' OF FADES IN RADIO TRANSMISSION 635 



REFERENCES 

1. Bullington, K., Inkster, W. J. and Durkee, A. L., Results of Propagation 

Tests at 505 mo and at 4090 me on Bevond-Horizon Paths, Proc. I.R.E., 
43, October, 1955, pp. 1306-1316. 

2. Rice, S. ()., Mathematical Analysis of Random Noise, B.S.T.J., 24, January, 

1945, pp. 46-156. 

3. Rice, 8. ()., Statistical Properties of a Sine Wave Plus Random Noise, B.S.T.J., 

27, January, 1948, pp. 109-157. 

4. Norton, K. A., Rice, P. L., Janes, II. B. and Barsis, A. P., The Rate of Fading 

in Propagation Through a Turbulent Atmosphere, Proc. I.R.E., 43, October, 
1955, pp. 1341-1353; Norton, K. A., Vogler, L. E., Mansfield, W. V. and 
Short, P. J., The Probability Distribution of the Amplitude of a Constant 
Vector Plus a Ravleigh-Dist'ributed Vector, Proc. I.R.E., 43, October, 1955, 
pp. 1354-1361. 

5. Palmer, D. S., Properties of Random Functions, Proc. Cambridge Phil. Soc, 

52, October, 1956, pp. 672-686. 

6. Favreau, R. R., Low, H and Pfeffer, I., Evaluation of Complex Statistical 

Functions by an Analogue Computer, presented at. the 1956 I.R.E. Con- 
vention. An abstract appears in Proc. I.R.E., 44, March, 1956, p. 388. 

7. McFadden, J. A., The Axis-Crossing Intervals of Random Functions, I.R.E. 

Trans, on Information Theory, IT -2, December, 1956, p. 146; also a paper 
scheduled for the March 1958 issue of Trans, on Information Theory. 

8. Pearson, Karl, ed., Tables for Statisticians and Biomelricians, Cambridge 

Univ. Press, 1931, Part II," Table VIII, Volumes of Normal Bivariate Sur- 
face, pp. 78-109. 

9. Kuznefsov, Stratonovich and Tikhonov, On the Duration of Exceedances of 

a Random Function, Journal of Technical Physics, 24, 1954. Translated 
from the Russian by Goodman, N. R., New York Univ. Coll. of Engineering, 
Scientific Paper No. 5, Engineering Statistics Group, March, 1956. 

10. Palm, Conny, Intcnsitalsschwankttngcn im Fernsprechverkehr , Ericsson Tech- 

nics 44, L. M. Ericsson, Stockholm, 1943. 

11. Ehrenfeld, S. and Goodman, R., private communication. 

12. Siegert, A. J. F., On the Roots of Markoffian Random Functions, Rand memo- 

randum, September 5, 1950. 

13. Wang, Ming Chen and Uhlenbeck, G. E., On the Theory of the Brownian 

Motion II, Revs. Modern Phys., 17, April— July, 1945, p. 323. 

14. Siegert, A. J. F., On the First Passage Time Probability Problem, Phys. 

Rev., 81, February 15, 1951, pp. 617-623. 



